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ABSTRACT 

We compute the number counts of clusters of galaxies, the logAMogS* relation, in 
several X-ray and submm bands on the basis of the Press-Schechter theory. We pay 
particular attention to a set of theoretical models which well reproduce the ROSAT 0.5- 
2 keV band \ogN-\ogS, and explore possibilities to further constrain the models from 
future observations with ASCA and/or at submm bands. The latter is closely related to 
the European PLANCK mission and the Japanese Large Millimeter and Submillimeter 
Array (LMSA) project. We exhibit that one can break the degeneracy in an acceptable 
parameter region on the Qq — o~$ plane by combining the ROSAT logiV-logS' and the 
submm number counts. Models which reproduce the ROSAT band logA^-logS 1 will 
have JV(> S) ~ (150 - 300)(SyiO~ 12 erg cm' 2 s" 1 )" 1 - 3 str" 1 at S > 10~ 12 erg cm" 2 s" 1 
in the ASCA 2-10 keV band, and A^(> S u ) ~ (10 2 - 10 4 )(^/100mJy)~ L5 str" 1 at 
S u £ lOOmJy in the submm (0.85mm) band. The amplitude of the logiV-logS' is 
very sensitive to the model parameters in the submm band. We also compute the 
redshift evolution of the cluster number counts and compare with that of the X-ray 
brightest Abell-type clusters. The results, although still preliminary, point to low 
density (Qq ~ 0.3) universes. The contribution of clusters to the X-ray and submm 
background radiations is shown to be insignificant in any model compatible with the 
ROSAT logAMogS. 

Subject headings: Cosmology — Galaxies : clusters of — Radio sources : extended - 
X-rays : sources 
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1. Introduction 

Clusters of galaxies are among the largest virialized structures in the universe and their 
importance as a cosmological probe is well-recognized. Thus their observations have been actively 
carried out in a variety of wavelengths including X-ray, optical, infrared, and radio bands. In the 
present paper, we focus on the theoretical predictions for the number counts of clusters of galaxies, 
the logiV-logS' relation, rather than more conventional statistics such as the X-ray temperature 
and luminosity functions (hereafter XTF and XLF) for several reasons; 1) temperature of X-ray 
clusters can be reliably determined only for luminous ones, and thus the statistics is inevitably 
limited, 2) such obtained XTF is to some extent weighted towards relatively rich clusters, and 
may be biased for the luminous species, 3) the XTF and XLF at high redshifts (z £ 0.1) are 
in fact model-dependent statistics, because the translation of the observed X-ray flux to the 
absolute luminosity, and of the observed number to the comoving number density can be done 
only by assuming specific values of the cosmological parameters (the density parameter, Qq, the 
dimensionless cosmological constant, Ao, and the Hubble constant Ho in units of 100 km/s/Mpc, 
h). 

On the other hand, the logAMogS relation is almost free from the above problems as long 
as the cluster identification (or separation from point-like sources) and the conversion of count 
rates to fluxes are reliable. Recent analysis of the ROSAT Deep Cluster Survey (RDCS, Rosati 
et al. 1995, 1997) and the ROSAT Brightest Cluster Sample (BCS, Ebeling et al. 1997a,b) 
has determined the logiV-logS' of clusters over almost four orders of magnitude in flux, i.e. 
5(0.5-2.0 keV) ~ 10~ 14 — 10~ 10 ergcm~ 2 s _1 . The number of identified clusters in the logiV-logS" is 
over 200, an order of magnitude larger than that for the commonly used XTF based on Henry <fe 
Arnaud (1991), and therefore the logAMogS data are statistically more reliable. 

Kitayama & Suto (1997, hereafter KS97) found that a set of cold dark matter (CDM) models 
reproduce the above R OS A T log iV-logS" data remarkably well over whole observed flux range, and 
simultaneously agree with the observed XTF and the COBE 4 year data. Nevertheless, there still 
exist some degeneracy of acceptable cosmological parameters. In the present paper, we explore 
possibilities to further constrain the models so as to break such degeneracy and discuss their 
implications in the following manner. 

First, we combine the logiV-logS" relations at different wavelengths, in X-ray and submm 
bands. The latter is of particular significance in relation to the future projects including the 
European PL A N CK mission and the Japanese Large Millimeter and Submillimeter Array (LMSA) 
project. The emissions from intracluster gas in X-ray and submm bands are originated from 
completely different physical mechanisms; the former is mainly due to thermal bremsstrahlung, 
and the latter is due to the inverse-Compton scattering of the cosmic microwave background 
(CMB) photons, i.e. the Sunyaev & Zel'dovich (1972, hereafter SZ) effect. As a result, the 
logAMogS' relations in these bands show very different parameter dependence. We note that 
the submm logA^-logS 1 was computed earlier by several authors (e.g., Barbosa et al. 1996; 
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Colafrancesco et al. 1997). Our analysis below differs from theirs in considering several CDM 
models consistent with the ROSAT log N-\ogS data, in including the relativistic correction to the 
SZ effect, and in making quantitative and extensive predictions for the number counts on the 
f^o — °8 plane. 

Second, we consider the cluster number counts incorporating redshift and/or temperature 
information in addition to the flux. In this way, we are able to discuss the evolution of cluster 
abundances on the basis of a cosmological model-independent and bias-free observable at high 
redshifts, which is in contrast to the approaches based on the XTF or XLF. We demonstrate a 
tentative comparison of our predictions with an observed sample from the ROSAT All Sky Survey, 
the X-ray brightest Abell-type clusters (Ebeling et al. 1996). 

Finally, we discuss the implications of our results for the X-ray background (XRB) and 
the submm background radiation (SBR). Since the logA^-logS" relation is closely related to 
the background radiation in the corresponding energy band, we may rigorously constrain the 
contribution of clusters of galaxies to the XRB and SBR. 



2. Number counts of clusters of galaxies in X-ray and submm bands 

2.1. X-ray flux from clusters 

Following KS97, we compute the number of clusters observed per unit solid angle with X-ray 
flux greater than S by 



oo 



N(> S)= dz d A (z)c 



o 



dt 

dz 



dS (1 + z) n M (M,z)-—— — - , (1) 

W J gas "-^band «>3 



where c is the speed of light, t is the cosmic time, d A is the angular diameter distance, T gas and 
-^band are respectively the gas temperature and the band-limited absolute luminosity of clusters, 
and hm(M, z)dM is the comoving number density of virialized clusters of mass M ~ M + dM at 
redshift z. 

Given the observed flux S in an X-ray energy band [E a ,E b ], the source luminosity L band at z 
in the corresponding band [E a {\ + z),Ef,(l + z)\ is written as 

L b and[K(l + z),E b (l + z)] = A7rd L 2 (z)S[E a , E b ], (2) 

where d L = (1 + z) 2 d A is the luminosity distance. Since the X-ray luminosity of clusters of 
galaxies depends sensitively on the details of the cluster gas density properties (distribution and 
dumpiness), its theoretical prediction as a function of M or T gas is difficult. In fact, a simple 
self-similar model predicts for the bolometric luminosity, L bol oc T g 2 as , which is inconsistent with 
the observed relation L bol oc T^ 3 ' 5 . Although a preheated cluster model might account for the 
latter (Kaiser 1991; Evrard & Henry 1991; Bower 1997), such theoretical models have not yet been 
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specified. Thus we adopt the observed L bol — T gas relation parameterized by 

L bol = L 44 (-^Y (1 + z)*: lO 44 ^" 2 erg sec- 1 . (3) 

As in KS97, we take L44 = 2.9, a = 3.4 and ( = as a fiducial set of parameters on the 
basis of recent observational indications (David et al. 1993; Ebeling et al. 1996; Ponman et al. 
1996; Mushotzky & Scharf 1997). Then we translate L hol (T gas ) into the band-limited luminosity 

-^band [^gas ; E\ , E<z\ &S 

^ba„ d [T gas , E a (l + z),E b {l + z)\ = L bol (T gas ) x /[T gas , E a (l + z),E b {l + z)], (4) 

where f[T gas , E\, E2] is the band correction factor which takes account of metal line emissions 
(Masai 1984) in addition to the thermal bremsstrahlung; the former makes significant contribution 
to the soft band luminosity especially at low temperature. Throughout this paper, we fix the 
abundance of intracluster gas as 0.3 times the solar value. Equations (||), @ and (||) relate S and 
T gas through Lband> and are used to compute equation ([[]). 

Assuming that the intracluster gas is isothermal, its temperature T gas is related to the total 
mass M by 

HrripGM 
7 3r vir (M, z f Y 

5.27(1 + z f ) ( V3 ( — ^ ) ^ nl /3 keV. (5) 

where k B is the Boltzmann constant, G is the gravitational constant, m p is the proton mass, // is 
the mean molecular weight (we adopt /i = 0.59), and 7 is a fudge factor of order unity which may 
be calibrated from hydro dynamical simulations or observations. The virial radius r v ; r (M, zA of a 
cluster of mass M virialized at Zf is computed from A v ; r , the ratio of the mean cluster density to 
the mean density of the universe at that epoch. We evaluate this quantity using the formulae for 
the spherical collapse model presented in Kitayama & Suto (1996b) and assuming for simplicity 
that Zf is equal to the epoch z at which the cluster is observed. 

Finally, we compute the mass function um(M, z)dM in equation (|l|) using the Press-Schechter 
theory (Press & Schechter 1974) assuming z = Zf as above. The effect of Zf ^ z is discussed 
by KS97 in this context, and the more general consideration of zj 7^ z is given in Lacey & Cole 
(1993), Sasaki (1994), and Kitayama & Suto (1996a,b). 

2.2. Submm flux from clusters due to the Sunyaev-Zel'dovich effect 

The inverse-Compton scattering of the CMB photons due to high temperature electron gas 
leads to distortion of the CMB spectrum (Sunyaev & Zel'dovich 1972). If the electrons are 



kp,T,r, 
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non-relativistic, the change in the CMB intensity observed at frequency v$ at a position angle 9 
from the cluster center is given by 



AI* R (x,9) 

y(0) 



2(fcgT CMB ) 3 -» 

————g{x)y{e), 
(hpc) 2 

x e x ( x 
^ (xcoth - -4 



(e 



30 k B T (z) .rf 
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oo Ttl e C 



(6) 
(7) 
(8) 



where x = hpu /(k B T CUBj Q), T CMBj0 = 2.726K is the present-day CMB temperature, hp is the 
Planck constant, cjt is the Thomson cross section, m e is the electron mass, and n e is the electron 
number density. At frequencies uq > 217GHz (or wavelengths < 1.4mm), g(x) becomes positive 
and galaxy clusters become positive sources. In what follows we consider mainly the submm band 
at the observed wavelength 0.85mm {vq = 350GHz) where the emission from clusters is fairly 
strong (x = 6.2 and g(x) = 6.7) and the ground observations are feasible. 

Then the total flux from a cluster located at redshift z is 



S™(x,M,z) 



Ai^ R ( x ,e)d 2 t 
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mJy, (9) 



where X is the hydrogen mass fraction for which we adopt X = 0.76 hereafter, and £lp is the 
baryon density parameter (we adopt Qp = 0.0125/i~ 2 ). It should be noted that the above flux 
corresponds to the value averaged over the entire cluster, and the observed logiV-logS' might be 
somewhat different depending on the details of the instruments specifically used in the survey. 
Aghanim et al. (1997) address this issue in detail. 

At T gas £ 10 keV, the electrons become relativistic and the above expressions need to be 
modified. We thus apply the relativistic correction to equation @ derived by Rephaeli (1995; 
see also Rephaeli & Yankovitch 1997). At 0.85mm, this leads to 4%, 11% and 16% reduction in 
flux at T = 3keV, 8keV and 12keV, respectively. We find that the relativistic correction at this 
wavelength is well fitted by 



AI, 



R 



1 - 0.013 



UkeVy 



AI, 



NR 



at 0.85mm. 



(10) 



While clusters with T gas ^ lOkeV are rare and the above correction does not make significant 
difference, we take it into account for completeness in computing submm luminosity functions and 
the logiV-logS . Since the above submm flux after the relativistic correction is explicitly expressed 
as a function of M and z, the submm logA^-logS 1 is computed in a similar manner to equation (|l]). 
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3. Number counts in X-ray and submm bands predicted in the cold dark matter 

models 

KS97 shows that CDM models with a certain range of parameters reproduce the ROSAT 
logiV-logS data as well as the XTF (Henry & Arnaud 1991) and the COBE 4 year data (Bunn & 
White 1997). To be definite, we hereafter consider five models with different sets of parameters 
summarized in table 1, and see if one can break the degeneracy of the models by comparing the 
logA^-logS relations in different bands. Throughout this paper, we assume that the primordial 
spectral index n is equal to unity and use the fitting formulae given in Kitayama & Suto (1996b) 
for the CDM mass fluctuation spectrum on the basis of Bardeen et al. (1986) transfer function. 

Figure |l] shows the logA-logS of clusters of galaxies in the soft X-ray (0.5-2 keV), hard 
X-ray (2-10 keV), and submm (0.85mm) bands. In the soft X-ray band, all the adopted models 
by definition reproduce well the observed logiV-logS' from the ROSAT Deep Cluster Survey 
(RDCS, Rosati et al. 1995, 1997) and the ROSAT Brightest Cluster Sample (BCS, Ebeling et 
al. 1997a, b). In the hard X-ray band, massive clusters make slightly larger contribution to the 
predicted logAMogS than in the soft band, and thus the lower Q,q models with flatter fluctuation 
spectra (i.e., more power at large scales) yield greater logAMog5 . All the models are shown to 
lie under the upper limit inferred from the number of all X-ray sources observed by ASCA at 2-10 
keV (Cagnoni et al. 1997), and they (except El) are also consistent with the number of clusters 
observed by HE AO 1 (Piccinotti et al. 1982) within the la errors. The predicted count in model 
El is smaller than the HE AO 1 result at the 3a level (possible incompleteness in the HE AO 1 
result would simply raise the observed data point and does not reconcile the discrepancy). 

In the submm band, on the other hand, the contribution from distant clusters becomes larger 
than in the X-ray bands (eqs and ||). Thus the lower models with slower evolution of 

fluctuations (i.e., more clusters at high z) yield greater logA-logS 1 as seen in figure [l|(c). Compared 
with the hard X-ray band, the degeneracy of the models are broken in much greater extent in 
the submm band. Therefore, future determination of the submm cluster counts will enable us to 
constrain the models more tightly. 

To see this clearly, we plot in figure || the contour maps of the cluster logiV-log5 in different 
bands on the f^o — °8 plane. For the 0.5-2 keV band, the 1 — a significance level derived from the x 2 
test between our theoretical prediction and the observations (see KS97 for detail) is plotted, while 
for the 2-10 keV and submm bands, the contours of the number of clusters per steradian (10 2 , 
10 3 , 10 4 from bottom to top) at 5(2-10 keV) = 10~ 13 erg cm 2 s x and 5^(0. 85mm) = 10 2 mJy are 
plotted respectively. The as values derived from the COBE 4 year data (Bunn & White 1997) is 
also shown for reference. Panels (a) and (b) show that the contours for the 2-10 keV band counts 
run almost parallel to the x 2 contour of the 0.5-2 keV band counts. In this sense, the future ASCA 
logiV-logS data will provide an independent consistency check of the ROSAT data. The shape of 
the submm logA-logS contours, on the other hand, is quite different, especially at high as, and 
thus should place complementary constraints on and a$. 
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In the conventional CDM power spectrum we have adopted, the spectral shape vary sensitively 
with CIq, h and Ob through the shape parameter T = Q,qH exp(— Ob — v / 2M^b/^o) (Sugiyama 
1995). In order to segregate the effects of changing the spectral shape, we further consider in the 
lower panels of figure Q the 'CDM-like' spectrum with the fixed shape parameter V = 0.25 inferred 
from the galaxy correlation function (e.g., Peacock 1996) independent of Oq, h and Ob- The 
corresponding logiV-logS' contours are similar to the conventional CDM case except at Oo ~ 0.2 
and Oo ^ 0.8, while the COBE normalized a% is very sensitive to the changes in the spectral shape. 

4. Evolution of cluster abundances 

The evolutionary behavior of the cluster abundance would provide a useful probe of the 
cosmological parameters as discussed by a number of authors (e.g., Viana & Liddle 1996; Eke, 
Cole, Frenk 1996; Oukbir & Blanchard 1997; Mathiesen & Evrard 1997). It is of particular 
importance because the current and near future observations will provide us with the rapidly 
increasing amount of information on high redshift clusters. Thus we consider the evolution of 
cluster distribution and discuss whether one can break the degeneracy among the models inferred 
from the soft X-ray logiV-logS' relation (table 1). 

Figure || exhibits the redshift evolution of the number of clusters in different bands. As 
expected, the evolutionary behavior strongly depends on the values of Qo an d erg', the fraction of 
low redshift clusters becomes larger for greater fio and smaller erg. It is indicated that one may 
be able to distinguish among these models merely by determining the redshifts of clusters up to 
z ~ 0.2 (see also discussion below). 

In order to characterize the physical properties of galaxy clusters, we have so far focused on 
their flux, because it is the simplest quantity which can be determined primarily from observations. 
Another major quantity of clusters is their temperature. Although the sample of clusters with 
measured temperature is still limited, an increasing amount of temperature information is going to 
become available from the ASCA (and future) observations. Furthermore, the temperature has an 
advantage in that it is insensitive to the detailed distribution of intracluster gas and much easier 
to model, compared to the X-ray flux or luminosity. This is in fact the main reason why many 
of the previous analysis of cluster abundance and its evolution are based upon the temperature 
rather than the X-ray luminosity (White, Efstathiou & Frenk 1993; Viana & Liddle 1996; Eke 
et al. 1996). However, the observed samples of clusters are usually flux limited, which needs to 
be kept in mind when one compares the theoretical predictions with the observed statistics. In 
our present framework, it is possible to incorporate explicitly the effects of limiting fluxes in the 
analysis based upon the cluster temperature. 

Figure || shows the cumulative distribution of clusters with temperature greater than T, 
redshift less than z, and with/without the X-ray flux limit. It is apparent that the presence of 
a flux limit affects significantly the temperature distribution and its evolution, especially at low 
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T and high z. If the flux limit is S = 10 -12 erg cm -2 s _1 in the 0.5-2 keV band, for instance, 
the number of clusters observed with T > 2keV and z < 0.3 per unit solid angle is an order of 
magnitude less than what would have been expected without any flux limit. 

The above discussion has clarified that the redshift evolution (at low z) and the temperature 
distribution of a given cluster sample with a specific flux limit can provide a useful probe of 
cosmological models. Unfortunately, we do not yet have a complete X-ray cluster sample with 
redshift and temperature information. However, one may still demonstrate the power of this 
approach using currently available 'best' sample. For this purpose, we adopt tentatively the X-ray 
brightest Abell-type clusters (XBACs, Ebeling et al. 1996), which is about 80 % complete and 
consists of 242 clusters with S > 5 x 10 -12 ergcm~ 2 s _1 in the 0.1-2.4 keV band and z < 0.2. 
We use the temperature and redshift data compiled in table 3 of Ebeling et al. (1996). About 
92 % and 30 % of all the clusters listed in this table have measured redshift and temperature, 
respectively. For the rest of the sample, the table also provides the redshifts estimated from the 
magnitude of the lOth-ranked cluster galaxy, and the temperatures estimated from the empirical 
L — T relation. Keeping in mind the incompleteness of the sample and uncertainties especially 
in the estimated temperature data, we simply intend to perform a crude comparison with our 
predictions. 

The results are plotted in figure H The sky coverage of the XBACs is hard to quantify 
mainly due to the uncertain volume incompleteness of the underlying optical catalogue as 
noted by Ebeling et al. (1996). In figure || therefore, we simply plot the real numbers of the 
XBACs and normalize all our model predictions to match the total number of the XBACs at 
its flux limit 5(0.1-2.4 keV) = 5 x 10~ 12 erg cm -2 s _1 and redshift limit z = 0.2. In general, the 
model predictions are shown to be capable of reproducing well the shape and amplitude of the 
observed distributions JV(> S,> T,< z) even when T and z are varied. Taking into account 
the incompleteness of the observed data and large statistical fluctuations at low numbers, the 
agreements with models L03 and L037 (both has Qq = 0.3) are rather remarkable. Since the 
shapes and amplitudes of the predicted curves in these figures are primarily determined by the 
value of Oq, this result provides a further indication for low universe. 



5. Contribution to the X-ray and submm background radiation 



The currently observed X-ray logiV-logS' especially at faint flux end (Rosati et al. 1997; 
Rosati & Delia Ceca 1997) places a model-independent constraint on the cluster contribution to 
the soft XRB. From the observed logiV-logS' relation in figure |l](a), the cluster contribution to the 
XRB is estimated roughly as 

dN 



I cl (0.5-2 keV) 
~ 2 x 10~ 10 I In 



S 



Sl2 



s n 



dS 
-I- 4.3 



dS 



erg cm 



~ 2 s x str 1 , 



(11) 
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where S12 = 10~ 12 erg cm -2 s" 1 , and we have fitted the 0.5-2 keV lc-giV-logS" by a broken power- 
law: N(> S) oc at S > S12 and iV(> S) oc 5~ L3 at S < S\2- The above estimate depends 
only weakly on the faint end flux 5 m ; n , and taking for instance S m i n = 10~ 15 erg cm -2 s _1 yields 
I c l (0.5-2 keV) ~ 2 x 10~ 9 ergcm -2 s _1 str _1 , which is less than 10 % of the observed XRB 
(Gendreau et al. 1995; Suto et al. 1996). 

To be more specific, we also compute numerically the contribution of clusters to the XRB. 
The XRB intensity I~-{Eq) from clusters at energy E$, as seen by an observer at present, is given 
by 



I (E ) = / dz y[ (12) 

vK 0J AttHJo (1 + z)y/n (l + zf - K(l + zf + A ' 1 ' 

where J u (E,z) is the comoving space-averaged volume emissivity at redshift z. We estimate this 
quantity by 

roc 

J u (E,z)= L v (E,M,z)n M (M,z)dM, (13) 
J 

where the luminosity L U (E, M, z) of a cluster of mass M at z is computed in a similar manner 



to Section |2.1| . It should be noted that the observed XRB spectrum is not derived from all 
sky surveys (in particular, at lower energy band), but from small regions in the sky where 
there are very few known bright X-ray sources. In order to compare with such observations, 
therefore, we also need to omit bright sources in our theoretical predictions by equation (|l^). 
For this purpose, we only consider contributions from clusters with flux below the critical value 
iS cr it(0.5-2 keV) = 10 -13 erg cm -2 s -1 , which roughly corresponds to the flux of the brightest X-ray 
source observed in the Lockman Hole by ROSAT (Hasinger et al. 1993). 

Figure ^a exhibits the XRB intensity from clusters of galaxies in the cases of models L03 
and El, both of which reproduce the observed logiV-logS' relation in the soft X-ray (0.5-2 keV) 
band. In these models, clusters contribute to the XRB only less than 20 % at Eq ^ 2keV, which 
is consistent with the rough estimate described above. Thus, clusters of galaxies cannot be the 
major sources for the XRB even in soft X-ray bands (Evrard & Henry 1991; Oukbir, Bartlett, 
Blanchard 1997). As is clear from the above analysis, this conclusion is based on the observed 
logiV-logS relation, and almost independent of the assumed cosmological parameters. In fact, if 
one adopts a theoretical L bol — T relation inferred from the self-similar assumption (Kaiser 1986) 
which yields a = 2 in equation (||), there exist some models in which clusters account for the entire 
soft XRB (Blanchard et al. 1992; Kitayama & Suto 1996a). Such models, however, are simply in 
conflict with both the observed L bol — T relation (a = 3 ~ 3.5) and the ROSAT logAMogS 1 . 

A similar analysis can be performed for the SBR. Since we do not yet have the observed 
log-/V-logS in the submm band, we fit the predicted logiV-logS of model L03 and obtain 

i/i?(0.85mm) ~ lO" 12 (^J~) ^ W m" 2 str~\ (14) 

where we have assumed a single power law with iV(> S u ) oc S^ 1 ' 5 . Taking S^mm = lmJy gives 
z/i^ (0.85mm) ~ 10" n W m~ 2 str -1 , which is about 3% of the detected SBR by Puget et al. 
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(1996). Results from the numerical integration shown in figure |6|b confirm that the clusters of 
galaxies contribute only less than ~ 5% of the SBR. 

6. Conclusions 

We have presented several cosmological implications of the number counts of clusters of 
galaxies. We have paid particular attention to the theoretical models which are in good agreement 
with the ROSAT logAMogS' in the soft X-ray (0.5-2 keV) band, and explored possibilities to 
further constrain the models from future observations in the ASCA hard X-ray and submm bands. 

In the submm band (0.85mm), models which reproduce the ROSAT \ogN—\ogS predict 
N(> S„) ~ (10 2 - 10 4 )(5 i ,/100mJy)- L5 str -1 at S u > lOOmJy. We have shown that the amplitude 
of the above relation depend sensitively on Qq and as, and in a substantially different manner from 
the ROSAT \ogN-logS . Thus, combining the two can break the degeneracy in the acceptable 
parameter region on the Qq — a$ plane. This indicates that the future observations by the 
European PLANCK mission and the Japanese LMSA project would provide powerful probes of 
these parameters. 

In the 2-10 keV band, the number counts show similar parameter dependence to those in 
the ROSAT 0.5-2 keV band, and we predict N(> S) ~ 200(5/10- 12 ergcm" 2 s -1 ) -1 - 3 str" 1 at 
S > 10~ 12 ergcm -2 s _1 in the ASCA 2-10 keV band. The ASCA logiV-logS' would therefore provide 
an important cross-check for our interpretation of the ROSAT log N-logS data. 

The evolutionary behavior of the number counts is also important to put additional 
cosmological constraints. We have exhibited that, given a complete flux limited cluster sample 
with redshift and/or temperature information, one can further constrain the cosmological models. 
We have performed a tentative comparison between our theoretical predictions and the recent 
compilation of the XBACs by Ebeling et al. (1996), which is the largest sample of galaxy clusters 
available to date. While the incompleteness of the sample and uncertainties in the temperature 
data still make it difficult to draw any definite conclusions from this comparison, it is interesting 
to note that our predictions reproduce well the evolutionary features of the XBACs and that the 
results, although preliminary, seem to favor low density (Qq ~ 0.3) universes. 

The cluster logAMogS also provides a tight constraint on their contribution to the background 
radiation in the corresponding energy band. Based on the logiV-logS' relation observed by ROSAT 
in the 0.5-2 keV band, we conclude that clusters of galaxies contribute at most ~ 20% of the total 
XRB and less than ~ 5% of the SBR. 
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Table 1. CDM model parameters from the ROSAT X-ray logAMogS" . 



Model 


n 


Ao 


h 




a 


7 


L03 


0.3 


0.7 


0.7 


1.04 


3.4 


1.2 


O045 


0.45 





0.7 


0.83 


3.4 


1.2 


El 


1.0 





0.5 


0.56 


3.4 


1.2 


L03 7 


0.3 


0.7 


0.7 


0.90 


3.4 


1.5 


L01q 


0.1 


0.9 


0.7 


1.47 


2.7 


1.2 
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Fig. 1. — The logAMogS relations of galaxy clusters for CDM models in (a) the soft X-ray (0.5- 
2.0 keV) band, (b) the hard X-ray (2-10 keV) band, and (c) the submm (0.85 mm) band. Lines 
represent the models listed in table 1; L03 (solid), O045 (dotted), El (short dashed), L037 (long 
dashed), and LOla (dot-dashed). Also plotted in panel (a) are the la error bars from the RDCS 
(Rosati et al. 1995, 1997), and the error box from the BCS (Ebeling et al. 1997a, b). In panel (b), 
the arrows indicate the logiV-logS' of all X-ray sources in the 2-10 keV band from ASCA (Cagnoni 
et al. 1997), and the error bar (1 a) is the number of clusters observed by HE AO 1 (Piccinotti et 
al. 1982). 
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Fig. 2. — Contour maps on the 0,Q-ag plane in (a) spatially flat (Ao = 1 — ^o) CDM models, (b) 
open (Ao = 0) CDM models, (c) spatially flat CDM-like models with the fixed shape parameter 
(r = 0.25), and (d) open CDM-like models with T = 0.25. In all cases, h = 0.7, a = 3.4, and 
7 = 1.2 are assumed. Shaded regions represent the la significance contours derived in KS97 
from the soft X-ray (0.5-2 keV) log./V-logS' . Dotted and solid lines indicate the contours of the 
number of clusters greater than S per steradian (10 2 , 10 3 , 10 4 from bottom to top) with S = 10 -13 
ergcm _2 s _1 in the hard X-ray (2-10 keV) band and with S u = 10 2 mJy in the submm (0.85 mm) 
band, respectively. Thick dashed lines represent the COBE 4 year result computed from the fitting 
formulae at 0.2 < ft < 1 by Bunn & White (1997). 
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Fig. 3. — Redshift evolution of the number of galaxy clusters. Upper panels show the cumulative 
number N(>S, <z) against z in (a) the soft X-ray (0.5-2 keV) band with S = 1CT 13 ergcm~ 2 s _1 , 
(b) the hard X-ray (2-10 keV) band with S = 10~ 13 ergcm~ 2 s _1 , and (c) the submm (0.85 mm) 
band with S u = 10 2 mJy Lower panels (d)-(f) are similar to (a)-(c) except for plotting the 
differential distribution dN(> S, < z)/dln z normalized by iV(> S). 
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Fig. 4. — Effects of limiting fluxes on the redshift and temperature distributions of galaxy clusters. 
Upper panels exhibit the number N(> S,>T,< z) against z for T = 2keV (thin lines) and T = 8keV 
(thick) in the cases of (a) S = (no flux limit), (b) S'(0.5-2keV) = 1CT 13 ergcm _2 s _1 , and (c) 
S'(0.5-2keV) = 10~ 12 erg cm -2 s -1 . Lower panels (d)— (f) are similar to (a)— (c) except for plotting 
N(> S, >T,<z) versus T in the cases of z = 1 (thin) and z = 0.1 (thick). 
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Fig. 5. — Tentative comparison with the XBACs (Ebeling et al. 1996). Upper panel (a) shows 
the log./V-log5 in the 0.1-2.4 keV band with the redshift limit of z = 0.2 (models L03 and O045 
almost overlap with L037 and El respectively). Middle panels exhibit iV(> S,>T,< z) versus z 
for 5(0. 1-2.4 keV) = 5 x 10~ 12 ergon -2 s-Mn the cases of (b) T = 2keV and (c) T = 5keV. Lower 
panels plot N(> 5, > T, < z) versus T for 5(0.1-2.4 keV) = 5 x 10~ 12 erg cm -2 s _1 in the cases of (d) 
z = 0.2 and (e) z = 0.1. All the model predictions are normalized to reproduce the total number 
of the XBACs at its flux limit 5(0.1-2.4 keV) = 5 x 10 -12 erg cm" 2 s _1 and redshift limit z = 0.2. 
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Fig. 6. — Contribution of clusters of galaxies to (a) the XRB, and (b) the SBR, in models L03 
(solid) and El (dashed). Also shown are the observed XRB data of ASCA (Gendreau et al. 1995), 
ROSAT (Hasinger 1992), ROSAT at 0.25 keV (Barber et al. 1996), and the tentative detection of 
the cosmic far- infrared background by COBE (Puget et al. 1996). 



